Driven-dissipative dynamics of a strongly interacting Rydberg gas 
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We study the non-equilibrium many-body dynamics of a cold gas of ground state alkali atoms 
weakly admixed by Rydberg states with laser light. On a timescale shorter than the lifetime of the 
dressed states, effective dipole-dipole or van der Waals interactions between atoms can lead to the 
formation of strongly correlated phases, such as atomic crystals. Using a semiclassical approach, 
we study the long-time dynamics where decoherence and dissipative processes due to spontaneous 
emission and blackbody radiation dominate, leading to heating and melting of atomic crystals as 
well as particle losses. These effects can be substantially mitigated by performing active laser cooling 
in the presence of atomic dressing. 
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I. INTRODUCTION 

Rydberg states are highly excited electronic states of 
atoms and molecules with large principle quantum num- 
bers n [IHS]. The remarkable properties of Rydberg 
states include their size r ~ a$n 2 with ao the Bohr ra- 
dius, implying huge polarizabilities a ~ n 7 and large elec- 
tric dipole moments between Rydberg states, and thus 
strong coupling to external electric DC and microwave 
AC fields. Rydberg states can be excited from atomic 
or molecular ground states by laser light via absorption 
of one or more photons, with Rabi frequencies scaling 
as ~ n -3 / 2 . On the other hand, Rydberg states are 
long-lived, with lifetimes r scaling as r ~ n 3 (r ~ n 5 ) 
for low (high) angular momentum states. For typical 
present experiments with n ~ 50 these can be as large 
as tens of /is. Furthermore, atoms prepared in Rydberg 
states interact strongly, and these interactions can be 
controlled and enhanced by external fields. In particu- 
lar, the van der Waals (vdW) interaction between re- 
states scales as VS/dw ~ (eao) 4 n 11 /r 6 , with e the elec- 
tron charge, and can be made both attractive and re- 
pulsive. In the presence of external DC or AC electric 
fields the Rydberg- Rydberg interaction is a dipolar inter- 
action, Vdd ~ (eaon 2 ) 2 /r 3 . The latter is both long-range 
and anisotropic. Additional tunability can be achieved 
using, e.g., Forster resonances. These phenomena have 
been explored in recent experiments [5U20|. 

The remarkable properties of Rydberg states, and the 
tunability and strength of interactions between Rydberg 
atoms is reflected in novel many-particle physics of Ry- 
dberg gases [2TT - f3T] . and provide the basis for applica- 
tions in quantum information processing, e.g. in im- 
plementing (fast) quantum gates [32 37 . An underly- 
ing principle is the Rydberg-blockade mechanism, as first 
proposed in Ref. [32] and [33]. In the blockade regime 
the presence of a single atom excited to a Rydberg state 



shifts the energy levels of the surrounding atoms within 
a characteristic radius rj, (~ /im), such that the laser ex- 
citation probability to the Rydberg state for any other 
atom within a volume ~ r 3 , is strongly suppressed. This 
excitation will be delocalized among all particles within 
ri, forming a "superatom" , with a collectively enhanced 
Rabi frequency y/Nt,£l. Here, N), is the number of atoms 
within r 3 constituting a superatom and f2 is the single 
particle Rabi frequency. Various phenomena in the fields 
of many-body physics, quantum information applications 
and quantum optics related to the dipole-blockade mech- 
anism have been discussed recently in Refs. [55144 7j . 

Present experiments on Rydberg gases as a many-body 
system created by laser excitation from BECs [JH US] , 
MOTs [6HT8 1 120 1 150]. atomic vapor cells [51] or optical 
lattices |19j have mainly explored the frozen gas regime. 
This corresponds to a short time dynamics, where the 
atomic motion can be ignored and (resonant) laser ex- 
citation of the Rydberg levels leads to large Rydberg- 
Rydberg interactions. In this limit atomic dissipation, 
e.g. spontaneous emission from Rydberg states, is negli- 
gible 52 and the dynamics maps to effective spin models, 
described by Hamiltonian dynamics [21] [23j |3J3 [53] . 

An alternative regime is the Rydberg- admixed gas 
regime [2"5"r[2"5l IM] : there the idea is to admix the strong 
Rydberg-Rydberg interactions by off-resonant laser light 
weakly to the atomic ground states, thus providing an 
effectively much smaller, but still tunable vdW or dipo- 
lar ground state dynamics. That is, instead of kDebye 
dipole moments of Rydberg states one obtains effective 
ground state dipoles in the range of tens of Debyes [55] . 
Thus kinetic energies and effective interactions can be- 
come comparable, while at the same time spontaneous 
emission from the Rydberg state due to off-resonant laser 
tuning is strongly reduced. This effectively extends the 
lifetime of the gas, making it reminiscent of the dynamics 
of dipolar gases of polar molecules, as reviewed in [5SH60j . 
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FIG. 1: (color online) Sketch of the proposed setup: Rydberg- 
dressed ground state atoms confined in the x-y-p\ane, e.g. by 
an optical lattice (not shown), are polarized perpendicular to 
the plane by a DC-electric field, Ed c , or an AC-microwave 
field, E mw (i). Atoms can either be in the dressed ground 
state \g) (green particles) with a dipole moment dg or in one 
of the intermediate Rydberg states \m) (red particles) with a 
dipole moment d m (see text), resulting in strong dipole-dipole 
interactions. 



An interesting question emerging from the above dis- 
cussion is, whether it is possible to form interesting con- 
densed matter phases with these engineered atomic dipo- 
lar gases; an example is provided by the formation of 
(stable) dipolar crystals of cold atomic gases in analogy to 
dipolar crystals discussed for polar molecules [61] . In con- 
trast to polar molecules [M for Rydberg-admixed gases 
dissipation and decoherence due to spontaneous emission 
is non-negligible for long times, and this raises the ques- 
tion of describing the long-time non-equilibrium dynam- 
ics of heating in such gases. 

Understanding the long-term dynamics of Rydberg 
dressed atoms can also have important applications be- 
yond many-body dynamics per se. For example, recent 
proposals have investigated the possibility to efficiently 
couple cold atomic Rydberg gases to comparatively hot 
molecular ensembles via long-range dipolar or vdW in- 
teractions in order to achieve Doppler and collisional 
Sisyphus-like |63j cooling schemes for molecules. 

The paper is organized as follows: In Sec.[ll]we give an 
overview of various atomic configurations studied, iden- 
tify the main questions, and summarize the main results 
of the paper. Technical details of our calculations can 
be found in the remainder of the paper: In Sec. 
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introduce the model and the notation which we will use 
throughout the paper. In Sec. |IV A| we discuss the Hamil- 
tonian for a single Rydberg-dressed atom in the presence 
of external static and electromagnetic fields. The in- 
teraction Hamiltonian of two Rydberg-dressed atoms in 
the presence of external fields is considered in Sec. |IV B| 
where Born-Oppenheimer potentials are derived for both 
DC-electric and AC-microwave fields. In Sec. IIV CI we 
study the validity of the 2D treatment. In Sec. [V] we 
derive Fokker-Planck equations for laser cooled ground 
state atoms in the presence of Rydberg dressing. We 
find an additional two-body diffusion term due to the in- 
teraction between two Rydberg-dressed atoms. The dy- 
namics of decay from the Rydberg state can be modeled 
by coupled Fokker-Planck equations for atoms effectively 




FIG. 2: (color online) Qualitative sketch of the energy levels 
(black lines), lasers (thin solid arrows) and decay paths (wig- 
gly arrows) for (a) the DC-electric field dressing scheme and 
(b) the AC-microwave dressing scheme. In both schemes the 
ground state \g) is weakly dressed with a Rydberg-state \r F ) 
(|r)) using a far detuned laser (red arrow) with Rabi frequency 
fi r and detuning A,- 3> Sl r . In the long-time limit cascaded 
decay from the Rydberg states \s) or \r F ) (gray wiggly ar- 
rows) will populate intermediate Rydberg states \m^ F )) due to 
spontaneous emission and blackbody radiation. In panel (a) a 
static electric field, Ed c = Fe z , polarizes the atoms leading to 
new Stark-split eigenstates, indicated with an index F. Both, 
the Rydberg state \r F ) and intermediate states \m F ) will ob- 
tain an intrinsic dipole moment d r and d m , respectively (blue 
thick arrows). In panel (b) a near resonant AC-microwave 
field, E mw (i), (yellow arrow) with Rabi frequency £l 3 cou- 
ples two Rydberg states \r) and \s) and induces an oscillating 
dipole moment proportional to d rs . The intermediate Ryd- 
berg states |m) have no dipole moment. 



in the dressed and laser cooled ground state and atoms 
in one of the intermediate excited states. In Sec. I VI Al we 
study numerically spontaneous emission and blackbody 
radiation of a single Rydberg dressed atom. Population 
of intermediate excited states after decay from the Ryd- 
berg state will lead to a fluctuating dipole moment and 
heating of the motion due to photon recoil. The case 
of an interacting ensemble of Rydberg dressed atoms is 
numerically studied in Sec. |VI C| We find that the fluc- 
tuations of the dipole moment lead to strong mechanical 
effects and heating of the external motion, which strongly 
depends on the atomic density and the dressing scheme 
used. We study numerically the effect of laser cooling 
and melting in the presence of an optical lattice with the 
aim to extending the lifetime of the gas. 



II. OVERVIEW OF RESULTS 

Before presenting a detailed discussion of the dynamics 
of strongly interacting Rydberg dressed atoms, we find it 
worthwhile to summarize the main features of Rydberg 
dressing in a gas of (realistic) atoms, and the resulting 
many-particle dynamics. The processes we summarize 
here are derived in detail in Sec. IIIII - IVII 

Motivated by proposals for many-body physics with 
dipolar gases [53 [M] , here we focus on the dynamics of a 
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Rydberg-dressed dipolar gas. We assume that the system 
is initially prepared in a crystalline state, e.g. by prepar- 
ing a Mott insulator state of atoms in an optical lat- 
tice, and turning on adiabatically the Rydberg admixture 
while switching off the optical lattice, so that a dipolar 
crystal is formed. Similar to the setups of self-assembled 
crystals for polar molecules studied in Refs. [5TJ [M] , we 
assume a two-dimensional (2D) configuration, i.e. that 
atoms are confined to the (x — y) -plane by an external op- 
tical field, e.g., an optical lattice, in the z-direction (see 
Fig. [I]) , such that in-plane dipole-dipole interactions are 
purely repulsive. As explained below, this will minimize 
collisional losses linked to the attractive part of dipole 
interactions. This analysis is readily extended to, and 
actually simplified in the case of isotropic repulsive vdW 
interactions in three dimensions [65 . 

The main question we want to address is how heating 
and dissipation affect the many-body dynamics due to 
decay from the Rydberg state in the long-time limit. In 
fact, for times long enough, spontaneous emission and 
blackbody radiation will inevitably redistribute the pop- 
ulation from the Rydberg state to various different ex- 
cited states, which invalidates the simple two-level ap- 
proximation for the internal dynamics of a single atom, 
mostly discussed so far [2TH3T] . Given the fact that atoms 
in these other excited states can interact very differently 
from those prepared in the dressed ground state, e.g. they 
can have different dipole moments, decay from the Ry- 
dberg state can lead to strong mechanical effects and 
collisional losses in an ensemble of interacting Rydberg 
dressed atoms. In the following we will investigate this 
complex many-body dynamics. In particular, we will 
explore how to mitigate and control these effects using 
laser-cooling or in-plane optical lattices for a realistic 
scenario where each atom comprises a large number of 
internal states. 



A. Atomic configuration 

The atomic configurations we have in mind are sum- 
marized in Fig. [2j (i) in panel (a) a DC-electric field 
with strength F polarizes each atom by splitting its en- 
ergy levels into the Stark structure; the new Stark-split 
Rydberg state \rp) obtains an intrinsic dipole moment d r 
which can be either parallel or antiparallel with respect to 
the external field, (ii) in panel (b) an AC-(microwave)- 
field of Rabi frequency f2 s is used to strongly couple two 
Rydberg states \r) and |s), which induces an oscillating 
dipole proportional to the transition dipole moment d rs . 

In both configurations of Fig. [2] we weakly admix the 
ground state \g) of each atom with the Rydberg state \rp) 
or |r), respectively, using an off- resonant continuous wave 
laser with Rabi frequency Q r and detuning A r 3> fi r . 
This immediately results in an effective dipole moment 
dg ~ (il r /2A r ) 2 do into the dressed ground state \g) ~ 
|p) + (Q r /2A r )\r), which can be tuned using the laser 

d r or do = d rs for 



DC-electric fields or AC-microwave fields, respectively, 
see Sec. 
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IV B This dressed ground state \g) has now a 



Deft comparatively long, lifetime ~ 1/Tg, where 



parameters, see Sec. IV Here, d { 



Tg ~ (fi r /2A r ) r r with r r the decay rate of the Rydberg 
state. 

For timescales which are comparable to, or even larger 
than 1/Tg population in each of the Rydberg states |rp) 
or |r) and |s) will be redistributed due to spontaneous 
emission and blackbody radiation among several excited 
states |m). In the long-time limit, after, e.g., a sponta- 
neous emission event from the Rydberg states, the atomic 
state will in general not return to the ground state di- 
rectly, but via a cascade process where several |m)-states 
are populated. Since each |m)-state has a finite lifetime, 
the cascade will not happen instantaneously. Population 
of these intermediate states can induce strong mechani- 
cal effects on the gas dynamics, leading to heating and 
losses, as explained below in Sec. |VI| It turns out that 
these effects depend crucially on how the dipole moment 
do in the Rydberg states is created, see Fig.ga) and (b). 

We analyze these effects in detail by performing semi- 
classical molecular dynamics simulations for an ensem- 
ble of interacting Rydbcrg-drcsscd atoms confined to a 
2D geometry including a large number of internal states. 
In Figs. [3] (a) and (b), we provide representative exam- 
ples for the cases (i) and (ii) above, respectively, for a 
system comprising N = 67 atoms in a box with hard 
walls. The external dynamics is treated classically, while 
quantum jumps from the Rydberg state account for spon- 
taneous emission and blackbody radiation, leading to a 
time-dependent dipole moment of the atoms. In the sim- 
ulation we use model-atoms with several thousands of in- 
ternal electronic states with the level structure and mass 
M of 85 Rb. At time t = the atoms are prepared in a tri- 
angular crystal structure with a given density n2D (see 
figure caption) at zero temperature T = 0. In Fig. |3j 
the thin dashed lines are the initial melting tempera- 
tures Tm = 0.089 d-n^p/kB as determined in Ref. [55] . 
with ks Boltzmann's constant. The solid blue and dot- 
ted green lines correspond to the mean atomic kinetic 
energy E^ n (t) and the total atom number N(t), respec- 
tively, plotted as a function of time t. In each plot, the 
simulation time corresponds to 1 /Tg , for the given choice 
of parameters. 



B. Rydberg dressing with a DC-electric field 

Figure [3]Ja) shows an example where we choose a very 
strong DC-field with strength F = 3 kV/cm in order to 
polarize 85 Rb atoms. Such a high electric field is moti- 
vated by studying mixtures of polar molecules and Ryd- 
berg atoms where both species are polarized by the same 
field [ST]. In order to avoid, e.g., field- ionization with 
threshold scaling as F mn ~ n~ 4 , the choice of Rydberg 
states is thus limited to those with a low principal quan- 
tum number, e.g. n ~ 16. Here, we use |r) = 1 16c?, m = 
0), with lifetime r r ~ 5.5 /is and field-induced dipole mo- 
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FIG. 3: (color online) Single trajectory of a semi-classical molecular dynamics simulation studying the non-equilibrium and 
melting dynamics of a 2D Rydberg-dressed crystal as a function of time t. Heating due to decay from the Rydberg states 
and population of intermediate states leads to a rapid increase of the mean kinetic energy, -Bkin(i) (solid blue line, left 

axis) and a reduction of the particle number due to losses (dotted green line, right axis). In panel (a) the dipole moment in 
the Rydberg state is created using the DC-electric field dressing scheme of Sec. |HB] and the same atomic parameters as given 



there. Initially we choose a density U2d = 1 /im 2 ■ Two insets show histograms of the momentum distribution at t = 27 /is and 



t — 40 fis. In panel (b) the dipole moment in the Rydberg state is created using the AC-microwave dressing scheme of Sec. II C 
with the same atomic parameters as given in the text. The initial density is n2D = 0.2 fim~ 2 . In both panels the simulation 
time corresponds to the effective lifetime of the dressed ground state atoms, r§ = 125 /is or Tg = 12.8 ms, respectively. The 
thin dashed line is the initial melting temperature, Tm, of the crystal. 



merit d r ~ 680 Debye. We choose a dressing laser with 
n r /A r = 0.42, result: 



125 fjs. 



ing m dg ~ 30 Debye and a lifetime 



The heating and loss processes described above are 
interaction-dependent (since the interaction strength 



Figure |3ja) shows that (i) (blue thick solid line) 
rapidly increases and exceeds Tm after about 10 /is. This 
is due to spontaneous emission events, and the ensuing 
population of various |m)-states, which have different 
dipole moments (see also scheme in Fig. [2ja)). Popu- 
lation of these intermediate states lead to a picture of 
fluctuating dipoles, where the energy of local crystal dis- 
tortions due to dipole-dipole interactions is rapidly redis- 
tributed among all particles as heat. The corresponding 
momentum distribution of the atoms is shown in the first 
inset. The heating effect from the fluctuating dipole mo- 
ments is dominant over single-particle recoil effects (with 
a characteristic rate of a few hundreds of kHz/ms, as 
discussed in Sec VI A below) . 

In addition, the blue line in panel (a) of Fig. [3] shows 
rapid large variations, or "spikes", of E^, e.g. at t ~ 
40 lis, which correspond to rare events, where one of the 
atoms populates a |m)-state with either a large dipole 
moment parallel to the DC-clcctric field and long life- 
time resulting in a strong repulsive interaction, or an 
anti-parallel dipole moment resulting in strong in-plane 
dipolar attraction. This out-of-equilibrium situation is 
reflected in a temporary change in the distribution of mo- 
menta of the atoms, as shown in the second inset by the 
appearance of small peaks at high-momenta. Such rare 
events result in collisional losses of the most energetic 
particles, similar to a filtering process, as well as small 
local "explosions" , where a large number of particles can 
be lost. 



scales as 



3/2 A 
'2D , 



and thus decrease considerably for a 



smaller atomic density. In Sec. |VIC| below we show how 
active laser cooling and an additional in-plane optical 
lattice can affect these processes. 



C. Rydberg dressing with a AC-microwave field 

In Fig. [3jb) we choose a strong and resonant AC- 
microwave field to couple the Rydberg states \r) = |50s) 
and |s) = \49p) (see level scheme of Fig. [2^b)), which 
have a transition dipole moment of d rs ~ 5.9 kDebye. A 
dressing laser with f2 r /A r = 0.14 yields dressed ground 
states atoms with dg ~ 30 Debye and Tg ~ 12.8 ms. 
The comparatively long lifetime is due to the choice of 
a state with a larger principal quantum number n ~ 50. 
Since F = intermediate Rydberg states |m) have essen- 
tially no dipole moment and we consider them as non- 
interacting. 

Figure ^b) shows that -EkinW increases with t. This 
is again due to fluctuations of the atomic dipole moment 
between the value dg and zero, corresponding to the atom 
being in the dressed ground state \g) or in one of the in- 
termediate states \m), respectively, together with photon 
recoil after a decay event. As discussed before, heating 
comes from the redistribution of the energy associated to 
local crystal distortions due to dipole fluctuations among 
all particles. For t < 1.5 ms the Rydberg dressed atoms 
are in a crystal-like phase, with E\^ n (t) < IzbTm- 
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In comparison with the DC-electric held case of 
Fig. [3^ a) the heating rate E^ n (T g )^ g of the AC- 
microwave-dressing scheme is approximately one order 
of magnitude smaller while the particle loss rate N(T g )j g 
is approximately equal. This is because the initial den- 
sity is smaller and the intermediate states have negligible 
dipole moments for AC-microwave scheme. 

In the remainder of this work we derive a microscopic 
model for the long-time heating dynamics of Rydberg 
dressed atoms and discuss quantitatively the dependence 
of the heating and particle loss rates on the system pa- 
rameters, e.g., the atomic density. We state under what 
conditions standard laser cooling can be performed in 
the presence of Rydberg dressing, and numerically inves- 
tigate how it counters heating effects. 



III. THE MODEL 

The purpose of this section is to introduce the phys- 
ical model and the notation that we will use through- 
out the paper. The system is illustrated in Fig. [lj N 
identical (alkali) atoms with momenta Pi at positions f j 
(1 < i < N) are confined to a quasi two-dimensional 
(2D) geometry in the (x-y)-plane by a strong confine- 
ment along the z-axis. 

The Hamiltonian dynamics of the system is governed 

by 



(ij) 
int i 



(1) 



i=l 



which consists of a sum over single-particle terms 
(see Sec. IV I and two-particle interaction terms 



H 



(ij) _ dj ■ d 3 - 3{dj ■ f)(dj ■ f) 



(2) 



which account for the dipole-dipole interaction between 
two atoms. Here, r = r^ — rj = f r is the relative distance 
between the atoms, the dipole operator of the i-th 
atom and eo is the vacuum permittivity. 



A. Internal level structure and setup 

We denote \a) — \n a ,£ a ,j a ,m a ) and hui a = E na i a j a 
as the unperturbed eigenfunctions and eigenenergies, re- 
spectively, of the atomic Hamiltonian H at . Here, n a is 
the principal quantum number, £ a the orbital angular 
momentum, j a the total angular momentum and m a pro- 
jection of the total angular momentum along a specified 
axis. 

In this manifold of states we focus on three specific 
states, see Fig. [4j (i) the energetic ground state of the 
atom denoted by \g) . All internal energies will be mea- 
sured relative to the energy of this state, e.g. hu> g = 0. 
(ii) a highly excited Rydberg states |r) with energy 



huj r and (iii) a lower- lying excited state |e) with energy 
huj e <C huj ri which will be utilized for laser cooling. All 
other states \a) can be occupied via spontaneous emis- 
sion or blackbody radiation from the high-lying Rydberg 
states. As we have already mentioned in Sec. [Hj the 
ground state of each atom is off-resonantly coupled to 
a high-lying Rydberg state \r) using a continuous wave 
laser with a Rabi frequency f2 r and detuning A r (^> f2 r ), 
see Fig. [3] Typically, in alkali atom experiments, the 
atoms are excited to the Rydberg state via two-photon 
transitions in which A r and £l r are the effective detun- 
ing and Rabi frequency. In addition, we use a laser held 
with Rabi frequency £l e and detuning A e to couple \g) to 
a nearby excited state |e) (Hui e <C Huj r ) forming a closed 
cycle for laser-cooling. 

Besides the general setup described above, we con- 
sider two different scenarios to create a strong dipole mo- 
ment in the Rydberg state: (i) an homogeneous electric 
DC field of strength F polarizes the atoms perpendic- 
ular to the 2D plane by splitting the energy levels into 
the Stark structure, see Fig. (4^a). The new Stark-split 
eigenstates \o.f) obtain a large intrinsic dipole moment, 
e.g. do = (ri?|d|rp). (ii) a resonant microwave held with 
Rabi frequency f2 s couples |r) to a nearby Rydberg state 
|s). Thereby the atom acquires a large oscillating dipole 
moment proportional to the transition dipole moment 
do = (r|d|s). In both scenarios, the dipole moments can 
be as large as kilo-Debyes. The corresponding dipole- 
dipole interaction between two atoms in the Rydberg 
states is governed by H lut of Eq. ^ . 



B. Master equation dynamics 

The dynamics of the driven-dissipative many-body sys- 
tem made of interacting Rydberg-dressed atoms is de- 
scribed by the following master equation 



p = -i [H, p] + Cp, 



(3) 



where the system density operator p acts both on the 
Hilbert space of internal \a\, a.2, . . . O/v) and external 
{(r 1; Pi), . . . (rjv, Pn)} degrees of freedom. This equa- 
tion is readily written down as an extension of the fa- 
miliar master equations of laser cooling for multi-level 
atoms including interaction between atoms [68l EH] . We 
denote L as the Lindblad operator, accounting for the 
decoherence due to spontaneous emission and blackbody 
radiation 



(4) 



i=l a, /3 



which results from the coupling between the atomic sys- 
tem and the vacuum modes of the electromagnetic field, 
which have been adiabatically eliminated. The decay rate 
from state |a) to |/3), 



r 



a/3 



p (SE) p(BBR) 
L a/3 T L a/3 ' 



(5) 
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is obtained as the sum of respective contributions from 
spontaneous emission (SE) and blackbody radiation 
(BBR) [70], with 



ri B fl BR) = n afj (T)A a (3, 



a/3 

.(SE) 

- a/3 



A a j3, U a > LOp] 
0, UJ a SC U)p. 



(6a) 
(6b) 



Here, A a p is the Einstein A-coefficient and n a p is the 
photon distribution, with 



3c 3 7rfe 



n a p(T) = 



1 



e H\u a -Ufi\/kBT _ 1 ' 

respectively. We use the notation 



(7a) 
(7b) 



d 2 kiV^(k)e- 



/3ot ' oc j3 



(8) 



to denote a general Lindblad term accounting for pop- 
ulation redistribution from level \a.i) to \Pi) of the i-th 

atom (1 < i < JV) with cr^ = |ai)(/3f|. With A Q(3 (k) we 
denote the angular distribution of spontaneous emission 
from level a to /3, which we assume to be a normalized 
and even function. The terms e ±li:a ' ik ' r ' describe the re- 
coil of the atom due to a spontaneously emitted photon, 
where k a p = (uj a — wp)/c is the wave number of the cor- 
responding transition, with c the speed of light and k a 
unit-vector in the direction of spontaneous emission. 



IV. HAMILTONIAN DYNAMICS 

We now turn to analyze the single particle and two- 
particle-interaction Hamiltonians of Eq. ([!]). In particu- 
lar, we derive the effective Born-Oppeheimer (BO) poten- 
tial surfaces for two interacting Rydberg-dressed atoms 
in the presence of either a DC-e lectric-h eld (Sec. IV B 1 ) 
or an AC-microwave-field (Sec. IV B 2). 



A. Single particle Hamiltonian 

The single particle Hamiltonian H^> of Eq. ([!]) consists 
of five terms 



H^{t) = |L + H$ + tf« + fl£>Jt) + H®,. (9) 

The first term in Eq. ^ accounts for the kinetic energy 
of the z-th atom with mass M. The second term 



(10) 
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FIG. 4: (color online) Energies uu a , states \a) (solid black 
lines), lasers (solid arrows) and decay paths (grew wiggly ar- 
rows) for (a) the DC-electric field and (b) the AC-microwave 
dressing scheme. Additionally to Fig. [2] a lower- lying excited 
state |e) is coupled in both schemes to the ground state \g) 
using a near-resonant laser (blue arrow) with Rabi frequency 
f2 e . State | e) decays directly to the ground state (blue wiggly 
arrow) with a rate T e yielding a closed cycle for laser-cooling. 
Cascaded decay from the Rydberg states r), \rp) (gray wiggly 
arrows) and \s) (not shown) will populate intermediate Ry- 
dberg states \m(F)) and |m'( F j) according to the decay rates 
T a p (see text). 



accounts for the internal atomic energy levels, in the ab- 
sence of external fields. The third term in Eq. ^ reads 



= • E dc = -d&F, 



(11) 



and describes the interaction of an atom with a static 
electric field, where dW is the atomic dipole operator 
of the i-th atom. The effect of a static electric field 
Edc = Fe z is to polarize the atoms along the field direc- 
tion e z , by splitting the energies into the Stark structure 
[71] . with F the strength of the field. The new Stark- 
split eigenstates \ap) have an intrinsic dipole moment, 
which in the linear Stark regime is approximately given 
by d a — (3/2)ea,Qn a (ni — n-i) with clq Bohr's radius and 
n\ and the parabolic quantum numbers of state \cxf)- 
In this regime the energy levels arc shifted proportional 
to the field strength, e.g. AE a = d a F. Coupling between 
adjacent n-manifolds can be neglected for field strength 



F < Fit, where -Fit 



is the Inglis- Teller- limit pQ. 
The term H^ eI in Eq. ^ describes the interaction of 
an atom with (e.g., microwave or optical) laser fields and 
consists of three terms 



ff£L(*) = ^ ; -w(*)+ ff . 



(i) 



a; -la 



(i) 
r-la 



where each term has the form 



H 



a/3 — laser 



(*) = 



Ml a 



00 



-i(k La ri — ui La t) 



(12) 



(13) 



Here k^ a = lo^/c is the wave- number of the laser with 
LJLa the frequency of the laser, h.c. denotes the her- 
mitian conjugate, and O a is the Rabi frequency. The 



first term H 



(i) 

eg— laser 



in Eq. (12) describes the coupling 
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of the atom to the cooling laser on the (|<?(F)H e (F)))- 
transition, with Rabi frequency Q e and frequency ojlc 
detuned by A e . The term H^_ laseI describes the cou- 
pling of the atom to the Rydberg-dressing laser on the 
(|(?(j7))-|r , (i?)))-transition, with Rabi frequency Q r and a 
frequency ujlv detuned by A r . In the following we are in- 
terested in the regime of large detuning A r 3> Sl r in order 
to weakly admix the Rydberg state to the ground state. 
Finally, i? sr _ laser describes the coupling to a microwave 
field strongly mixing the Rydberg-states |r) and |s), with 
Rabi frequency Q s and laser frequency u>ls detuned by 
A.. 

In addition H^ p of Eq. ^ accounts for external trap- 
ping potentials. 



B. Two-particle Hamiltonian: Born-Oppenheimer 
potentials 

In this section we study the interaction between two 
Rydberg dressed atoms in the presence of an external 
(Sec. " 



static electric field 
(Sec 



IV B 1 ) or a microwave field 



IV B 2 ) . For relative distances between the atoms 



a Q nf, 



larger than the size of the Rydberg atom r 
and in the presence of external fields the atoms interact 
via dipole-dipole interaction governed by the Hamilto- 
nian of Eq. Q. In particular, we derive the BO poten- 
tial surfaces which, in the adiabatic approximation, play 
the role of effective interaction potentials [6T| (Eqs. ( 18 ) 
(28)). They form the basis for the analysis of the time- 



dependent dynamics of Rydberg-dressed atoms which we 
discuss below in Sec. IVIIandlVl 



1. DC-electric field 

We consider two atoms in the presence of a static DC- 
field of strength F, driven by the Rydberg-dressing laser, 
in the configuration of Fig. |4]jb). The microwave and 
the cooling laser are absent, i.e. f2 e = £l s = 0. The 
new Stark-split eigenstates \g F ) ~ \g) and \r F ) are ob- 
tained by diagonalizing iJ at + i?dc and the detuning of 
the dressing laser A r is defined relative to the shifted 
energy levels. To obtain the BO potentials we first ne- 
glect dissipation and treat the position operators f , as 
parameters r,. Within this limit each atom can be de- 
scribed by a two-state model consisting of \g) and |t\f), 
coupled by a dressing laser. This is valid for distances 
larger than r n ~ (D/AEn) 1 / 3 , where diabatic crossings 
between BO-surfaces of neighboring n-manifolds can be 
neglected. Here, D = c?Q/(47reo) is the dipolar coupling 
strength and AE n ^> hA r is the energy separation be- 
tween neighboring states. In a rotating frame the sin- 
gle particle Hamiltonian describing this model system 
reduces to 



H® = -hA r a 



(i) 



\ y rq 



(14) 



where the operator cr£l = \a F ) (/3f\ acts on the new 
Stark-split eigenstates, i <E {1,2}, and A r = — oj r 
is the detuning from the ( \r F ) — \g) )-resonance. Position 
dependent phases of Eq. |l3| ) will be included in Sec. [V 
where they lead to recoil kicks from laser absorption anc 
spontaneous emission, and Doppler shifts. 

Since the DC-electric field aligns the dipoles of the 
atoms along the direction of the field and the dominant 
dipole moment is do = (rp|d|rF), the interaction term 
in Eq, (|| is 



(y) _ D(l - 3 cos 2 i?) 



(15) 



where ?? is the angle between the dipole axis and the 
radial vector between two atoms. When the atoms 
are confined in a 2D-plane, e.g. by a strong optical 
field, the angle is fixed to i? = tt/2 resulting in a 
purely repulsive interaction. Within this model, the 
total Hamiltonian [Eq. ([IJ] for two atoms in the basis 
{\g,g) , \ r F,g) , \g,r F ) , \r F> r F )} reads as 



/ 



h = n 



-A r 



\n r o 
o 



2 o 

-A r 







\tt r \n r V{r)-2A r 



(16) 



with 



hV{r) = {r Fl r F \H^ ] \r F ,r F ) = 



J 



D 



(17) 



The BO potentials for the scattering of two Rydberg- 
dressed atoms are obtained by diagonalizing the Hamil- 
tonian of Eq. (161 for fixed relative position and zero 
kinetic energy, see Fig. [5] We assume that the linewidths 
of the corresponding states are smaller than the energy 
separation between them, and that the relevant kinetic 
energies are small enough, such that Landau-Zener tran- 
sitions between different BO surfaces can be neglected 
(secular approximation). In this case the resulting po- 
sition dependent eigenvalues act as potentials in each 
state manifold [STJ [75]. Asymptotically (e.g., at large 
distances where the dipole interaction is negligible), the 
new dressed eigenstates are \g) — |p) + (fi r /2A r ) |t>) and 
\f F ) = \r F ) — (Sl r /2A r ) \g). We note that the dynam- 



ics governed by the Hamiltonian of Eq. ( 16 ) for atoms 
initially in the ground state is restricted to the symmet- 
ric subspace made of the three states which asymptoti- 
cally connect to the states \g,g), (\f F ,g) + \g,r F ))/V2, 
and Iff, fir), with energies slightly perturbed by the 
(Rydberg-dressing) laser field; the antisymmetric state 
(\r F ,g) — \g,r F ))/i/2 is decoupled, and does not con- 
tribute to the dynamics. 

In the following, we focus on blue detuning, e.g. A r > 
0. Fig. [5] shows the BO-potcntials as a function of 
the interparticle distance r for a specific choice of pa- 
rameters. Because of the choice of blue-detuning, the 
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FIG. 5: (color online) Energy eigenvalues E(r) (dressed BO 
potential surfaces) of Rydberg-dressed atoms confined in a 
2D geometry obtained by diagonalizing the Hamiltonian of 
Eq. (161. Here, r = r(cos <p, sin if, 0) is the 2D coordinate in 
the plane with z — 0. Atoms are polarized by the DC-electric 
field dressing scheme of Sec. |IV B 1| and A r = 2Q r . Energy 
surfaces are labeled using arrows: Egg (solid green), E g ?+? g 
(dashed-dotted) and E,-f (dashed). At the Condon point r = 
r c an avoided crossing leads to a rapid change of the ground 
state interaction potential E gg : for r > r c atoms prepared in 
the dressed ground state \g) are weakly interacting, E gg ~ 
(n r /2A r ) 4 do/r 3 , while for r < r c the potential inherits the 
character of the Rydberg-Rydberg interaction, E gg ~ do/r :i . 



figure shows that the ground state BO-potential Egg 
(solid green line) corresponding to the energy of the two- 
particle asymptotically in the state \g, g) has the highest 
energy. The energy of this BO-potential can be calcu- 
lated perturbatively up to fourth order in the small pa- 
rameter Q r -c \V(r) — 2A r \ as 



E gg (r) 



2A r 

^_ 
2A r 



4A3 
8A3 I 2A 



4A2 (y - 

; / r, x 4 



2A r ) 
V(r). 



(18a) 
(18b) 



where the second line is valid in the limit V(r) <C A r , 
i.e. r 3> r c (see below). While the first two terms 
on the r.h.s. of Eq. (18) are simple light shifts, the 



equation shows that at large distances the effective 
ground state BO-potential has an effective spatial 
dependence V gg (r) — (fl r /A r ) 4 d r ~/r 3 . As explained in 
Refs. [25 28, 54 , this means that by dressing the parti- 
cles with the laser field we have achieved a dipole-dipole 
interaction for atoms prepared in their (dressed) ground 
state, and the strength of the interaction is tunable by 
varying the ratio f2 r / A r . 

Figure [5] shows that the energy of the dressed two- 
particle ground state state E gg (solid green line) is 
strongly affected by the dipole-dipole interactions, and 
for V(r) ~ 2A avoided crossings occur among the BO- 
potentials of all symmetric states. Note that the energy 
of the antisymmetric state (\tf9f) ~ Iqf^f))/" 1 /^ is un- 
coupled and not shown. In particular, there is a resonant 



4 



87ren/iA T 



1/3 



(19) 



between the ground state BO-potential and other energy 
surfaces. As a consequence, there is a sudden change 
in the slope of the energy surface for r ~ r c , where 
the ground state BO-potential inherits the character 
of the one that asymptotically connects to the energy 
of the state \fF,fp), and becomes strongly repulsive. 
This effect has been discussed in Refs. [61] [65l [72] 
in the context of so-called blue-shielding techniques, 
where the strong repulsion for r < r c does not allow 
for particles to come close to each other in a scattering 
event, thus preventing collisional losses due to, e.g., 
collision-induced ionization at short distance. In the 
reminder of this work, we will be mostly interested in 
confining the dynamics to distances r > r c . 



2. AC-microwave-field 

In this section we consider the Hamiltonian dynam- 
ics of two atoms in the presence of a linearly polarized, 
near-resonant microwave field with Rabi frequency Cl s 
coupling the Rydberg states \r) and \s), see Fig. Qb). 
The ground state is again weakly dressed with the state 
\r) using an off-resonant laser with Rabi frequency Q r 
and a large detuning A r 3> O r . The DC-electric field and 
the cooling laser are absent, i.e. 51 e = and F = 0. The 
explicit choice of a near-resonant microwave field is to ob- 
tain large dipoles for atoms in the dressed ground state, 
which scale as ~ (fl r /A r ) 2 , similar to the DC-electric 
field case of Sec. |IVB 1[ Again neglecting dissipation 
and the external degrees of freedom for a moment we are 
left with a three level system consisting of the ground 
state \g) and the two Rydberg states \r) and \s). Since 
the electric DC-field is absent, these states arc the bare 
eigenstates of H at . Below we show that, by a judicious 
choice of systems parameters, it is possible to obtain an 
interaction strength and Condon radius of similar mag- 



nitude as in the previous scheme of Sec. IV B 1 above. 

In a frame rotating with the laser frequencies the single 
particle Hamiltonian of Eq. ([lj for this model system 
reduces to 



hA r o 
hVL 
~2~ 



(0 



- h(A r 



(<) 

rg 



gr , 



~2 



(20) 



where the operator o^L = |a) (/3\ acts on the bare eigen- 
states, i € {1,2}, and A s = a;^ s — lu s is the detuning of 
the microwave laser from the (|r) — |s))-resonance. The 
corresponding dipole-dipole interaction Hamiltonian of 
Eq. ([2| reduces to 

3 cos 2 ?? r 



to) 



47T60 



rs sr 



(21) 



9 




Is) 
■1+) 



with the corresponding eigenenergies 



£4 



-A, 



(24) 



FIG. 6: (color online) Illustration of the AC-microwave dress- 
ing scheme: (left) energies uj a , states \a) (solid black lines) 
and lasers (solid arrows) of the system described by the Hamil- 
tonian of Eq. ( |20[ ) . The unitary transformation U rs of Eq. ( 22 1 
diagonalizes the Rydberg subspace {\r), \ s)}. (right) In a ro- 
tating frame the new states |+) and |— ) with energies E± 
are separated by an energy difference V^l + A| and individ- 
ually coupled to the ground state with Rabi frequencies fi± 
according to the Hamiltonian of Eq. ( 25 1 . 



After the transformation, the single particle Hamilto- 
nian of Eq. (20) in the basis {\g) , |+) , |— )} is 




(25) 



with d defined as before. 

In the case of a near-resonant microwave field, where 
A s <C Q s it is convenient to perform a unitary transfor- 
mation [73] 

C/W = exp I \ tan- 1 (g) («# - <#) } , (22) 

which diagonalizes the Hamiltonian in the subspace 
{\ r )i ) l s )il- The corresponding new eigenstates are 
l±)i = a T l s )i ±a ± |r)i with 

a± = 4= I 1 ± / As I , (23) 



with the effective Rabi frequencies Q± — a±Q r . In this 
transformed picture the new eigenstates |±) are coupled 
to the ground state by lasers with Rabi-frequencies il± 
and have a dipole moment d± = (±|d|±) = ±(r|d|s) = 
±d , see Fig. [6] 

Again, the dynamics of the symmetric and 
antisymmetric states are decoupled. For two 
particles the dynamics of the symmetric sub- 
space is governed by the Hamiltonian H = 
+ + Hj„ t which, represented in the basis 

{|fl,fl>,^(lfl,+> + l+,5»,^(lfl,-> + |-,fl»,^(|- -) + 
l+, +», -) + !-> +)), 4KK -) - l+, +))}, reads 



{ 
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-2A r 





-n. 





V - 2A r y 



(26) 



where we assumed exact resonance of the microwave field, 
i.e. A s = 0. 



The BO-potentials are obtained by diagonalizing the 
Hamiltonian H sym , leading to the new dressed eigenstates 
which parametrically depend on r. An analytic expres- 
sion for the BO-potential of two atoms in the dressed 
ground state can be obtained perturb at ively in the limit 



n r < {A r ,n s } 



P MW 



2A r C% 
4A? - «2 



1 - 



n 2 r (AAl + 3Ql 
(4A2 - W s f 



- 16A^ - (4A 2 . + 3ftg)(»g - 4A 2 .)] 

(4A2 - ^2)2 ( 2 y Ar _ 4A 2 + fl 2) 



(27) 



Note, that there are two resonances at A r = ft s /2 and 
2l/(r c MW )A r - 4A 2 . + n 2 s = 0. The first one corresponds 
to the level crossing between the states \g—) and | ) 
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FIG. 7: (color online) (a) Energy eigenvalues E(r) (dressed 
BO potential surfaces) of Rydberg-dressed atoms confined in 
a 2D geometry obtained by diagonalizing the Hamiltonian of 
Eq. ( |26[ |. Here, r = r(cos ip, sin tp, 0) is the 2D coordinate 
in the plane with z = 0. Atoms are polarized by the AC- 
microwave dressing scheme of Sec. |IV BT] and A r = 4fi r , fl s = 
1.5A r and A s = 0. Energy surfaces are labeled using arrows. 
The energies 

and £™5(r) are strongly affected by 
the dipole interaction and at the Condon point r^f w = g(a)r c 
an avoided crossing leads to a rapid change of the ground state 
interaction potential £^ w (r) (thick green line): for r > r c 
atoms prepared in the dressed ground state \g) are weakly 



interacting, Eg 



f(a)(n r /2A r ydt,/r- i , while for r < r c 



the potential inherits the character of the Rydberg-Rydberg 
MW - dl/r 3 . The figure shows that the Condon 



interaction, E gg 



radius is slightly shifted to larger values compared to the DC- 
dressing scheme of Fig. [5] Panel (b) shows the dependence 
of the interaction strength, /(a), (green and solid line) and 
the Condon radius, g(a) (grey and dashed line) defined in 
Eq. ( 28 1 on a. There exists a region Aq where /(«) > 1 and 
g(a) < 1.5. 



(red detuning) or \g+) and \-\ — h) (blue detuning) with 
the two-particle ground state \gg). The second one cor- 
responds to a Condon point r^ w , similar to the one we 
have discussed before for the DC-electric field case. 

We now focus on the case O s /2 < A r in which the 
single-particle |+)-state with energy E + ~ — A r + ~Q S ps 
— e (0 < e <C A r ) gets almost degenerates with the 
ground state \g), while the state |— ) with energy E- ~ 
—A r — |fi a s» — 2A r is separated by a large energy gap of 
2A r . In order to investigate the behavior of Eq. ( 27 ) near 
this resonance we set 



The latter inequality comes from the fact that we assume 
il r is the smallest frequency scale for the non-degenerate 
perturbation theory to be valid. Therefore, these two 
states must not be exactly degenerate. With this as- 
sumption we find for the ground state energy and the 
Condon radius 



=E +f(a)[^) Vir, 



.MW 



"o 



(28) 



87reo^A r ' 



with £ const . = 0?/[2A r (l - a 2 )] - nj(l + a 2 )/[8A3(l - 
a 2 ) 3 ], f(a) = 2a 2 /(l - a 2 ) 4 and g(a) = 1/(1 - a 2 ) 1 / 3 . 
In the limit a — > the state \s) is not coupled to \r) and 
we obtain the same light shifts as in Eq. (l8| ). Note that 
in this limit the interaction strength vanishes, f(a) — > 0. 
This is due to the fact that in the absence of a DC-electric 
field the bare state \r) has no intrinsic dipole moment. 

Figure^a) shows the BO-potentials for the symmetric 
states as a function of the interparticle distance r for a 
specific set of parameters near resonance, as discussed 
above. The energy of the dressed two-particle states 

| ) and |++) is strongly shifted by the dipole-dipole 

interaction, and avoided crossings occur among the BO- 
potentials of all symmetric states. This leads to a sudden 
change in the slope of the energy surface of, e.g. E gg , at 
the Condon radius r^ w . 

The functions f(a) and g{a) versus a are shown in 
Fig. [7^h). Increasing a towards 1 will on one hand in- 
crease the effective ground state interaction potential ac- 
cording to f{a). On the other hand it will increase the 
Condon radius r^ w according to g(a). Fig. fab) shows 
that there is a region, Act, for which f(a)> 1 but 
1 < g(a) < 1.5 between 0.45 < a < 0.84. Operating 
in this region leads to formal similar interaction strength 
and Condon radii as for the DC-electric field dressing 
scheme, e.g. Eq. (18) and Eq. (19 1. 



C. Validity of the 2D treatment 

In this section we examine in detail under what criteria 
we can treat the system of interacting Rydberg dressed 
atoms as purely 2D in nature. As mentioned before, the 
atoms are trapped in the (x — y)-plane by a strong har- 
monic confinement along the z-direction. The resulting 
three dimensional (3D) potential in relative coordinates 
reads 



V ZD {v)= l -Mujlz 2 +hE(v 



(29) 



where the first term is the harmonic confinement with 
trapping- frequency ui±. The second term, E gg (r), is the 
BO potential of two interacting dressed ground state 
atoms, obtained by numerically diagonalizing the Hamil- 
tonian of (16), taking into account the full 3D character- 



a2A r with a < 1 — n r /A r . istic of the dipole-dipole interactions. 
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FIG. 8: (color online) Contour plots of the 3D-potentials 
V 3D (r) and V^'lr) of Eq. ([29} and (30} are shown in panels 
(a) and (b), respectively, in units of Vb = 2A r (f2 r /2A r ) 4 . 
Here, r = (p cos ip,p sin <p,z) is the relative distance between 
two atoms, with p the in-plane radial coordinate and z the 
transversal coordinate. Brighter regions represent stronger 
repulsive interactions. Two saddle points (red circles) located 
at (/J*, ±2*) separate the repulsive from the attractive short- 
range region. 



Fig. [8la) shows a contour plot of the 3D poten- 
tial V3d[t) of Eq. (291 in the (p-z)-plane where p = 
\/ x 1 + y 2 . We consider 85 Rb atoms with A r = 2ir x 
250 MHz, fl r = 2vr x 100 MHz, r c = 520 nm and 
uj± = 2ir x 200 kHz. In the figure, darker color corre- 
sponds to deeper potentials. The potential exhibits two 
saddle-points, at (p*, z*) = (0.78, ±0.39) r c , with a height 
Vsd(p*, £*) = 220 pK ■ ks, which serve as an energy bar- 



rier separating the repulsive long-range dipole-dipole in- 
teraction regime from an attractive short-distance regime 
[6T1 [64] . For relative kinetic energies smaller than the 
height of the potential barrier, the interaction is purely 
repulsive and the system can be stabilized against col- 
lapse due to the attractive part of the interaction. For 
blue detuning of the dressing laser there is a resonant 
Condon point at p = r c for z = 0, where we observe a 
rapid increase of the interaction potential, which is dis- 
cussed in Sec. 



IVB1 



Along the axial direction z (for 
p = 0) there is no resonant point and the potential ap- 
proaches smoothly zero where inelastic or reactive colli- 
sions will occur. 

For relative distances larger than the Condon radius, 



r c [Eq. (19[)], the interaction can be approximated by 
dipole-dipole interaction. In this case, the 3D potential 
can be rewritten as 



KZ 2 



p 2 - 2z 2 



7» 



(30) 



which is reminiscent of polar molecules [5TJ [M]. Here, 
z = z/r c and p = (x 2 + y 2 ) 1 ^ 2 jr c . The corresponding 
contour plot is shown in Fig. [8]^b) . The dimensionless pa- 
rameter k characterizes the strength of the confinement 
relative to the interaction 



2A r 



Mj\rl 
8hA r 



r c 

a ho 



(31) 



where a^ — WfijMuJx is th e harmonic oscillator length 
in transversal direction. For k » 1 the confinement in 
z-direction dominates over the dipole-dipole interaction. 



Analytic analysis of Eq. ( 30 ) yields the saddle points 



±22* 



and 



/5k 1 / 5 ' 



(32) 



The height of the potential at the saddle point, corre- 
sponding to the height of the energy barrier separating 
the repulsive region from the attractive region, is 



V, 



^\z„P,) 



V 



3/5 



(33) 



which is independent of A r for a fixed ratio of f2 r /A r . 
For the same parameters used above we obtain k = 30 
which yields a potential barrier at the saddle point cor- 
responding to = 150 pK. Comparing Fig.[8]Ja) and (b) 
one can see the rapid change of the interaction strength 
at the Condon radius for Eq. (29 1 (panel (a)) leads to 



a repulsive shield which leads to a slightly higher po- 
tential barrier corresponding to = 220 pK compared 
to the case of a pure inverse-cubic interaction potential 
(panel (b)) with T* = 150 pK. Additionally we find that 
the saddle-point in panel (a) is shifted to larger p values 
than in panel (b). 

In a 3D scenario two particles undergo collapse when 
they have large enough collisional energy to overcome the 
energy barrier at the saddle points. In the numerical sim- 
ulations of Sec. |VI| we will treat the system as purely 2D 
with a l/r 3 -potential (see Eq. (42 1), assuming particles 
are being lost if they collapse in 3D. For a specific set 
of parameters we use the latter analysis to calculate the 
height of the potential at the saddle point, V3d(p*, 
in a 3D scenario. The height of the saddle point will then 
be translated into a critical 2D distance, pi oss , at which 
Egg(p ioss ) = V 3D (p±, z±), with E g g defined in Eq. (fl8])(b). 



V. LASER COOLING OF DRESSED RYDBERG 
ATOMS 

In this section we examine in detail the dissipative pro- 
cesses in Eq. ([3} and study laser cooling of interacting 



12 



atoms in the presence of Rydberg dressing. As an ex- 
ample we consider Doppler cooling but the model can be 
extended to sub-Doppler cooling schemes. For simplic- 
ity of the analytic treatment in Sec. |V A| we will first 
consider model atoms where we neglect decay from the 
Rydberg state in order to derive semiclassical Fokker- 
Planck equations which describe the cooling dynamics in 
the presence of Rydberg-dressing and interactions. In 



Sec. |VB| we will include decay from the Rydberg state 
which couples the dressed ground state via rate equations 
to one of the intermediate states. Again for simplicity of 
the analytic treatment in this section, we will first con- 
sider model atoms with a single intermediate state |m). 
In Sec. |VI| we will numerically investigate the interplay 
of laser-cooling and heating due to population of inter- 
mediate states of an ensemble of dressed Rydberg atoms 
taking into a large number of internal states. We show 
that laser cooling can alleviate the heating dynamics for 
interacting Rydberg atoms described above, extending 
the lifetime of strongly interacting phases in these sys- 
tems, as for self-assembled crystals. 



A. Fokker-Planck equation 



In the following we will derive equations of motion 
for the external dynamics (position and momentum) of 
model atoms consisting of a ground state \g) which is 
coupled with a far detuned laser to a Rydberg state |r) 
with Rabi frequency fl r and detuning A r 3> £l r - Addi- 
tionally, the ground state \g) is coupled to a lower-lying 
excited state |e) using a counter-propagating laser with 
Rabi frequency fl e and detuning A e . The state |e) de- 
cays directly to the ground state \g) with a fast decay 
rate r e , realizing a closed cycle transition. Atoms which 
are both in the Rydberg state interact via dipole-dipole 
interaction described by the Hamiltonian Eq. (15 1. 



The derivation of equations of motion for the external 
degrees of freedom is done in four steps [69, 74J: (i) first 
we use the Wigner function formalism to map the density 
operator p onto a quasi-probability distribution in phase 
space 
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(34) 



(ii) We derive the equation of motion for W using the 
master equation Eq. ^ and expand the resulting equa- 
tion of motion up to second order in terms of the pho- 
ton momentum in order to obtain a positive probability 
function in the semiclassical limit, (iii) We adiabatically 
eliminate the Rydberg state \r) in the limit of a large de- 
tuning A r . The resulting two-level systems - consisting 
of a dressed ground state \g) — \g) + (Q r /2A r )\r) with 
dipole moment dg and an excited state |e) - interact only 
when the atoms are in the dressed ground state, (iv) Fi- 
nally, we adiabatically eliminate the fast internal degrees 
of freedom of this effective two-level system in favor of 
the (much slower) external dynamics. Due to the Rabi 
oscillations between the dressed ground state \g) (with 
dipole moment dg) and the lower-lying excited state |e) 
(with negligible dipole moment) the dipole of the atom is 
fluctuating in time. This leads to an additional diffusion 
term in the equation of motion in addition to the stan- 
dard Doppler-cooling diffusion terms. After a lengthy 
calculation one obtains for two atoms in their dressed 

( 2") 

ground state, W~~' = (gg\ W\gg), the following equation 



of motion 




which is valid when A r 3> 0, r ,il e ,T e 3> hk^ g /m, where 
k eg is the wave number of the (|e)-|g))-transition. Addi- 
tional to the free motion term on the left hand side, the 
first term on the right hand side describes single particle 
Doppler cooling with the cooling force f c = /3ki e kx e • p 
and the standard diffusion operator 1>i = D^&^ k ac- 
counting for spontaneous emission and diffusion due to 
the cooling laser in various spatial directions k € {x, y, z} 
|74j . Here we assumed two counter-propagating laser 
beams and expanded the resulting radiation pressure 
forces up to second order in ki e . The third term of 



Eq. ( 35 ) proportional to 



dv 

4A2(F-2A r ) 2 <9n 



dEgg 



(36) 



accounts for the interaction between the two atoms. We 
note that this term can be rewritten, using the potential 



of the dressed ground state, Eq. ( 18 ), obtained in section 
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IV B Besides those standard laser cooling terms, [59ll74j . 
the term 



4^r 



(r* + 4A») 



fl2 



_d d_ 

dpi <9p 2 



(37) 



is a two-body diffusion term which accounts for the fluc- 
tuations of the force (which is present only in the dressed 
ground state), due to Rabi-oscillations between |e) and 
\g). For near-resonant laser light it is proportional to the 
population and lifetime of the excited state. In the limit 
fi2 <C fir e k£ e this diffusion term is small compared to 
the single-particle diffusion terms and one can approxi- 
mate the system with two interacting atom, which are 
independently laser cooled. 



B. Quantum jumps and rate equations 

A quantum jump occurs when the atom in the Ryd- 
berg state \r) decays to an intermediate state \m) due 



to spontaneous emission or blackbody radiation. In the 
following, we allow for the possibility that the interme- 
diate state |m) is long-lived, with a lifetime ~ 1/T mg 
comparable with the external dynamics and thus cannot 
be adiabatically eliminated. This will lead to two cou- 
pled equations of motion: (i) an equation for laser-cooled 
atoms in the dressed ground state, with dipole moment 
dg = (il r /2A r ) 2 d r , and (ii) an equation of motion for 
atoms in the \m) -state (and possibly with a dipole mo- 
ment d m ) 7 which do not experience a cooling force. 



To simplify the notation for the analytic treatment, 
we consider only a single atom in the presence of the 
cooling and dressing lasers. The generalization to two 
atoms is straightforward. After preforming steps (i)-(iv) 
of the latter section we find the following coupled FPEs 
for the dynamics of the dressed ground state \g) and the 
intermediate state Ito) 



dt M dr / 9 







d 



r ^ ■ — ) 
dt M dr " 



-TgW ( s 1] + (T m +V mg )w2 ) + 



r w {1) 

1 TO m 



r 



2A 7 



dp 



■ f c + £>i 



ir. 



(i) 



w. 



(i) 



(38a) 
(38b) 



where w£ ] = (a\W^\ a) is the single particle Wigner 
function for the atom in state \a) (a £ {g,m}). The 
equations are coupled by the effective ground state decay- 
rate 



2A r 



(39) 



due to optical pumping via the Rydberg state \r). The 
terms T> a p are standard diffusion operators |74) account- 
ing for spontaneous emission from state |a) to |/3). 



VI. DISSIPATIVE DYNAMICS OF DRESSED 
RYDBERG ATOMS 

We will now analyze numerically the dissipative pro- 
cesses of Eq. Q , which are associated with the finite life- 



time of excited Rydberg states. In Sec. VIA below, we 
consider the case of a single atom, in which we investigate 
the population of intermediate states due to spontaneous 
emission or blackbody radiation. As a consequence, the 
atom acquires a time-dependent dipole moment. In addi- 
tion, spontaneous emission and blackbody radiation act 
as small heating sources due to the photon recoil. In 



the driven-dissipative dynamics of laser-cooled and in- 
teracting Rydberg-dressed atoms as a prescription for a 
molecular dynamics simulation. The effect of population 
in intermediate states will be more substantial in the case 
of an ensemble of Rydberg-dressed atoms, as we discuss 
in Sec. 



VIC 



In particular, we find that the dominant 
heating effect originates from the time-dependent dipole 
moment which induces strong mechanical effects in the 
many-body system. A crucial point is that, the time- 
dependence of the dipole moment exhibits characteristi- 
cally different behavior in the DC-electric field case com- 
pared to the AC-microwave one. For example, in the case 
of a static electric field (F > 0) the states lap) can have a 
large parallel or antiparallel dipoles, d Q = (ai?|d|aF), on 
the order of kilo-Debye, causing strong dipole-dipole in- 
teractions with other Rydberg-dressed atoms. Whereas, 
in the case of a AC-microwave field (and no DC-electric 
field) the intermediate states do not have any dipole mo- 
ment. 



A. Decoherence of a single atom 



The dressed ground state \g) has a finite lifetime 



Sec. VI B we will summarize our previous discussion of r g = T~ , where T g — (fl r /2A r ) T r , and l/r r is the 
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FIG. 9: (color online) Relative decay rates r r , n ^/r r of the 
|r) = |16d 5 /2, mj = ~5/2) state of 85 Rb. Green and gray 
bars indicate decay rates to \np)- and |n/)-states, respectively. 
In both cases we summed over j and rrij. The upper and 
lower panel show the relative decay rates for T — K and 
T — 300 K, respectively. The total decay rate of Tied = 2nx 
51 kHz (43 kHz) for T — 300 K (T = 0) agrees well with 
calculations carried out in [70]. Note that in the upper and 
lower panel the y-axes are cut at 5 % and 10 %, respectively. 
For states with a higher decay rate percentage numbers on 
the left side of the bars indicate their value. 



lifetime of the Rydberg state. The latter is in general 
given by blackbody radiation and spontaneous emission, 
which redistribute population from the \r) (|rF))-state to 
all possible intermediate states \m). Such a decay event 
from the Rydberg state is followed by a cascade process 
where several intermediate states |m) can be populated. 
It is crucial to note that, the cascade process in general 
does not happen instantaneously, due to the finite life- 
time of the intermediate states. As we discuss below this 
has far reaching consequences in the long-time dynamics 
of Rydberg-dressed crystals. 

We calculated the decay constants T a p given in Eqs. 



( 6a ) and ( 6b I numerically using both quantum defect the- 
ory [751177] and a model potential method [78]. With 
the first method we calculated the decay matrix in- 
cluding all angular momentum states up to n = 18, 
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FIG. 10: (color online) Stark map: Atomic energy levels 
E a (F) of 85 Rb as a function of the field strength F around 
the state 1 16d, m = 0) . Note that we plot only states with 
magnetic quantum number m = 0. States with an angu- 
lar quantum number £ > 3 have a negligible quantum de- 
fect and show a linear Stark effect, while s, p and d states 
show a quadratic. The energy of the state which for F — > 
connects to 1 16d, m — 0) (thick green line) is well separated 
in energy from neighboring states up to a field strength of 
F w 300 kV/m. 



while with the second method we obtained the decay 
matrix up to n = 110 including s, p, d, / and g- 
angular momentum states. With T a = "Y^a T a p we de- 
note the total decay rate of the state \a) and To = Tg = 
(Q r /2A r ) 2 T r = (tt r /2A r ) 2 J2p F r /3 is the effective decay 
rate of the dressed ground state. Fig.|9]shows the branch- 
ing ratio of decay in the absence of external fields from 
the 1 16(2) state to \np) and \nf) states, where we summed 
over j and rrij levels. For T = K (upper panel) the only 
contribution to the decay rate comes from spontaneous 
emission, which favors decay to low-lying states. About 
60 % of the population in the 1 16<i) state decays to the 
\5p) state. The lower panel of Fig. [9] shows the branching 
ratio for T — 300 K. In this case the decay rate is de- 
termined not only by spontaneous emission but also by 
blackbody radiation, which in addition leads to substan- 
tial decay to neighboring states. Hence, only 51 % of the 
population in the |16d) states decays directly to the \5p) 
state. From this it is clear that in current cold atom ex- 
periments performed in room temperature environments, 
the role of high-lying intermediate states cannot be ne- 
glected in the long-time limit. 

The internal dynamics corresponding to a cascade pro- 
cess of the electron towards the ground state following a 
decay event vim the Rydberg states is given by 
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(40) 



where p a is the probability of being in the state |a(_F)) 
in the absence (presence) of a static electric field, F = 
(F > 0). Note that in the case of a DC-electric field 
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FIG. 11: (color online) Electric dipole moment d a = 
(ot f | d z \ ocf} of the Stark split states |qf) of a 85 Rb atom in 
an electric DC-field with strength F — 3 kV/cm. States are 
numbered with increasing energy, e.g. the state 1 16d, m — 0) 
corresponds to a — 1017. The figure shows that depending 
on the state the dipole moment d a can be positive (parallel) 
or negative (antiparallel) leading to a repulsive or attractive 
interaction. 



(F > 0) the decay rates r QF( g F are a sum of several decay 
rates between bare states, T a p, according to the contri- 
butions of various different angular momentum states to 
the Stark split states \ap) and \/3p). 

DC-electric field: The Stark structure of atomic states 
(including only m = states) is shown in Fig. 10 for a 
electric field along e z . States with an angular momentum 
I > 3 are almost degenerate and exhibit a linear Stark 
effect, while lower angular momentum states, which have 
a larger quantum defect which breaks the degeneracy, 
exhibit a quadratic Stark shift. Fig. [10| shows that the 
\16d) state (thick green line) is well separated from the 
neighboring states by an energy gap of ~ 200 GHz up to a 
field strength of F ~ 3 kV/cm and experiences a strong 
energy shift, corresponding to a large dipole moment. 
For an electric field E<j c = F e z only m = states are 
coupled and the resulting dipole moments are polarized 
along the direction of the external field. Fig. [TT] shows 
the z-component of the dipole moment, d Q , for different 
states |op) and a field strength of F = 3 kV/cm. The 
figure shows that states \ap) gain a large dipole moment 
of hundreds of Dcbye which can be positive and negative. 
The dipole moment of the state which for F — > connects 
to \16d,m = 0) is about 680 Debye. 

The time evolution of Eq. ( 40 ) is readily simulated. We 
prepare the atom in the dressed ground state \g), that is 
Pg — Po = 1- After a time step At much smaller than 



the timescales associated with the decay rates a random 
number determines if and to which internal state 
the atom makes a transition (quantum jump) according 
to the rates of Eq. (40 1 . During each time step the dipole 



moment d(t) of the atom is equal to the dipole moment 
of the current state, rf Q(F) . This new configuration then 
again propagated another time step At. 

Figure [T2|a) shows an example result for the time 
evolution of the atomic dipole moment d(t) of a 85 Rb 
atoms polarized by a DC-electric field with strength 




FIG. 12: (color online) Time evolution of the atomic dipole 
moment d(t) according to population of various internal states 



governed by Eq. (401 for a Rb atom. In panel (a) the atom 



is polarized using the DC-electric field dressing scheme while 
in panel (b) it is polarized with the AC-microwave dressing 
scheme using the same parameters as given in Sec. |VI A| In 
both cases the dipole moment of the dressed ground state is 
dg — 30 Debye. In the case of the DC-electric field dress- 
ing scheme (panel (a)) the dipole moment fluctuates between 
large positive and negative values due to population of inter- 
mediate states, while in the case of the AC-microwave dressing 
scheme (panel (b)) the dipole moment jumps between d(t) = 
and d{t) = 30 Debye. 



F = 3 kV/cm. The ground state is coupled to the 
Rydberg state \rp) which in the absence of a DC-electric 
field (F — 0) connects to the |r) = \16d, m = 0)-state 
using a Rydberg-laser with (Q r /2A r ) — 0.21. We 
initially prepare the atom in the dressed ground state 
\gp), and compute the time evolution using Eq. (40 1. 



The red trajectory in Fig. [I2[a) shows a typical result 
for the deviation of the atomic dipole from the dipole 
moment of the dressed ground state d g — 30 Debye 
as a function of time t, in units of the ground state 
lifetime r g = 125 fis. The figure shows large positive 
and negative spikes for the values of d, followed by long 
times where the system is in the dressed ground state, 
with dipole moment dg. The large "spikes" correspond 
to the population of intermediate states \cxf)i during 
the cascade process towards the ground state, following 
a spontaneous emission event. These large positive and 
negative fluctuations of d(t) will cause strong dipole- 
dipole interactions and hence large mechanical effects 
when a gas of interacting Rydberg atoms is considered. 

Microwave dressing and F=0: In this case we cou- 
ple the ground state 1 5s) of Rubidium to the state 1 17s) 
using a laser with (n, r /2A r ) 2 = 0.05. An additional 
microwave field couples the states 1 17s) and \17p) such 
that the dressed ground state obtains a dipole moment 
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* g — 30 Debye. The blue trajectory in Fig 
shows a typical trajectory for the time evolution of 
the dipole moment as a function of time in units of 
Tg = 90 /is. For this dressing scheme the intermediate 
states \m) do not possess a dipole, and thus the dipole 
fluctuates between long periods when it has the value d g , 
corresponding to the atom in the dressed ground state, 
to periods where the atom has no dipole, corresponding 
to the cascade processes following spontaneous emission, 
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FIG. 13: (color online) Single kinetic energy trajectories, -Ekin, 
as a function of time t for a single Rydberg-dressed 85 Rb atom. 
We observe an increase in kinetic energy due to photon recoil 
due to decay events from the Rydberg state. The thick green 
line is a linear fit of 50 single trajectories (thin black lines) 
which gives a mean heating rate of 2-zr x 106.3 kHz/ms. 



like a "blinking dipole". 

In comparison, in the DC-clcctric field case (blue tra- 
jectory in panel (a)) the dipole fluctuations are much 
larger and can take both positive and negative values. 
Whereas in the microwave case (red trajectory in panel 
(b)) the dipoles fluctuate between zero and dg. 

In the following we examine the external dynamics of 
a single atom. It is dominated by small momentum fluc- 
tuations associated with the photon recoil following a se- 
ries of cascade decay events. As an example, we simulate 
heating due to decay from the \r) = \ l6d) state using the 
same parameters as before. The result of a molecular dy- 
namics simulation is shown in Fig. |13[ where we plot sin- 
gle kinetic energy trajectories as a function of time (thin 
black lines) . We propagate the internal dynamics accord- 
ing to Eq. (40). Each decay event is associated with a 



momentum kick in a random direction corresponding to 
the transition, see next Section. It is shown in the figure 
that as time progresses the kinetic energy increases due 
to photon recoils. We perform a statistical average over 
50 runs of simulations (thick green line) which yields an 
average heating rate of 2ir x 106.3 kHz/ms. It will be 
shown below that for high enough densities this single 
particle heating rate is much smaller than the heating 
rate due to interactions associated with fluctuations of 
the dipole moment in the many-body discussed 
in the next section. 



B. Molecular dynamics simulation 

In our semiclassical description the state of each atom 
at time t is specified by Pi(i), a,}, where r^i) and 

Pi(t) are the center-of-mass coordinate and momentum 
of the i-th atom, respectively, and a € {<?, m, m', . . .} de- 
notes the internal electronic state with i — 1, . . . , N. The 
internal state is either the ground state \g), dressed by 
the admixure of the Rydberg state, and the excited state 



\e) from the laser cooling (see Fig.|4|, or one of the many 
intermediate states \m) populated during a decay cascade 
back to the ground state. The dynamics of the gas is de- 
scribed by a set of coupled equations, as written down for 
the case of one or two atoms in Eqs. ( 38 ) and ( 35 ) , respec- 



tively. These results can be generalized immediately to N 
atoms. Mathematically, they are a set of coupled Fokker- 
Planck equations for the external motion, and rate equa- 
tions for the internal electronic states, describing a com- 
bined diffusive and jump Markov process for the prob- 
ability density W^ N '(ri, p 1? a\\ . . . ,r N ,p N , a^'yt). This 
stochastic process is readily simulated. 

Consider atom i, which we assume to be in the dressed 
ground state \g). Its center of mass motion obeys the 
Langevin equations 



Pi 

M ' 

Ef(ji) 
1 int 



/3 Pl + f W + pW 



(41a) 
(41b) 



where on the right hand side of Eq. (41b) we sum over 



the forces from all the other ground state atoms, as well 
as the forces from atoms in one of the intermediate states 
\m) , with 



d ai d aj d 1 
47re dvi |r< - Tj 



(42) 



Here, d ai = dg is the effective dipole moment of the 
dressed ground state (see Eqs. (18) and (28) for the DC- 
and AC-field cases, respectivelyfand d a . can either be 
equal to dg if the j-th atom is in the dressed ground state 
or equal to d m if the j-th atom is in the state |m). In 
the DC-electric field case d m = (mp\d\mF) , while in the 
AC-microwave field case d m = 0, see Sec. VIA In addi- 
tion, we have added in Eq. ( |41b[ ) the familiar terms de- 
scribing possible laser cooling [74|. and an external trap- 
ping force fix • The last term in Eq. (41b) is a stochastic 
force 79 from quantum fluctuations due to the recoil 
of spontaneous emission events from both the Rydberg 
state repopulating the ground state, but also from laser 
cooling, obeying 



(F^(t)F^(t'))=D k S kl d ij 5(t-t'), 



?<J) 



(43) 



with k, I £ {x, y, z} and i,j € 1, . . . , N and D k the dif- 



fusion coefficients of Eq. ( 35 ) . We make the simplifying 



assumption that the cross-noise term of Eq. (37) aris- 



ing from the fluctuating ground state dipole due to laser 
cooling discussed in Sec. [V]is negligible. 

Similarly, atom i in one of the intermediate states \m) 
obeys the equation of motion 



Pi 
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1 int I 



(44) 



where again f^' 
Note that atoms in one of the (intermediate) Rydberg 



are the forces of Eq. ( 42 ) and d a 
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states |m) are assumed to be neither trapped 
nor laser cooled, /3 = 0. 

Optical pumping from \g) to one of the intermediate 
states | to) and the following cascaded decay, \g) — > \m) — > 
. . . — > \m!) — > \g), back to the dressed ground state will 
redistribute the atomic populations according to the rate 
equations (40). We can simulate this many-body dynam- 
ics by starting with a given atomic configuration, and 
propagating Eqs. (41 1 and (44) for a time step At, much 



smaller than the timescales corresponding to the above 
mentioned rates. We then determine the probabilities for 
atom i in state CKj to make a transition (jump) to another 
internal state according to the rate equations ( 40 ) for op- 



tical pumping and decay, and pick a new configuration 
according to these probabilities. Each cascade of quan- 
tum jump is associated with momentum kicks Kk rm , 
ftk m >g in random spatial directions according to the dis- 
tribution N a p of Eq. p| and with k a) g the wavevector 
of the corresponding transition. This new configuration 
is then again propagated another time step according to 
equations Eqs. (41) and (44). 



Due to population of intermediate states in the long- 
time limit and the resulting time-dependent dipole mo- 
ments we expect a rapid heating caused by fluctuations of 
the interaction force described above. In addition we also 
expect atomic losses. For example, when a ground state 
atom collides with an atom in an intermediate | a) -state 
with a large dipole moment the kinetic energy released 
in the collision can (i) be of the order of the height of 
the saddle point (see Sec. IV C). In this case the atoms 



will overcome the potentia barrier and experience an at- 
tractive interaction. Or, (ii) the energy released can be 
of the order of the depth of external confinement in the 
transversal direction. A third loss mechanism is be Pen- 
ning ionization at distances smaller than 4n 2 ao (5UHS2- 
Each of these processes leads to two-body or one-body 
losses. In the simulations, we utilize a model including 
two-particle losses, where two particles are lost whenever 
their relative distance r becomes smaller than a certain 
critical radius r\ oss (see Sec. IV C). 



C. Effects of decoherence on the dynamics of many 
interacting Rydberg atoms 

In this section we numerically investigate the effects 
of spontaneous emission and blackbody radiation on the 
long-time dynamics of an ensemble of interacting Ryd- 
berg atoms by performing semi-classical molecular dy- 
namics simulations. As discussed above population of 
intermediate states |to) with different interaction prop- 
erties, e.g. different dipole moments, will lead to strong 
fluctuations of the interparticle forces in the many-body 
case. 

Using semi-classical molecular dynamics simulations 
we study the resulting non-equilibrium dynamics of an 
ensemble constitution of N — 67 Rydberg-dressed 85 Rb 
atoms. Initially they are prepared in a perfect triangu- 
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FIG. 14: (color online) DC-electric field case: We numerically 
investigate the heating rate (upper panel) and the remaining 
particle number after a time t§ = 125 fis (lower panel) as a 
function of the density U2d for different laser cooling param- 
eters /3 and lattice depths Vq. Parameters see text. 



lar crystal in a box with hard walls at T = 0. For both 
dressing schemes we choose system parameters such that 
the ground state dipole moment is dg — 30 Debye, corre- 
sponding to a melting temperature of Tm — 0.6 /iK (2ttx 
12 kHz) for U2D = 1 /^m~ 2 . Additionally, the effect of 
laser-cooling and a triangular in-plane optical lattice is 
analyzed in order to stabilize the atomic crystal. Note 
that, we assume only dressed ground state atoms experi- 
ence both the cooling laser and the lattice with a depth of 
50 En, with Eji — 2tt x 3.8 kHz, while atoms in one of the 



intermediate states \m) do not, see Eqs. (41) and (44) 



1. DC-electric field dressing 

Fig. [14] analyzes heating in a gas of Rydberg dressed 
atoms using the DC-electric field scheme of Fig.|4^b) with 
the same parameters as in Sec. |II B| For A r = 2n x 
250 MHz the Condon point is at r c ~ 520 nm. For a 
transversal trapping frequency of u± = 2ir x 200 kHz 
the height of the potential barrier at the saddle point is 



150 /iK. Each marker in Fig. 14 is an average over 50 
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FIG. 15: (color online) AC-microwave field case: We numer- 
ically investigate the heating rate (upper panel) and the re- 
maining particle number after a time r§ = 12.8 ms (lower 
panel) as a function of the density U2d for different laser cool- 
ing parameters /3 and lattice depths Vq. Parameters see text. 



runs of the simulation, where we simulated the dynamics 
for a time t„. 



Fig. 14 'a) shows the heating rate, jgE^ n (t = r§) as a 



function of the atomic density for different cooling and 
lattice parameter, with -E"kin(^) the mean kinetic energy 
of all N atoms. In the case of no laser cooling, /3 = 0, 
there is a strong dependence of the heating rate on the 
density, while for a cooling rate of (3 = 2tt x 30 kHz the 
heating rate is almost zero and a steady state is realized 
on a timescale /3 _1 <C Tg. This steady state is due to the 
interplay of laser cooling and heating together with loss 
of high energy particles from the 2D confinement. This 
is shown in Fig. |14[ b), where the number of remaining 
particles after a time Tg is shown as a function of the 
atomic density. Again, there is a strong dependence on 
the density and also on the laser cooling rate j3. For 
densities n 2 D ^ 0-04 /im" 2 the effect of the fluctuating 
dipoles can be neglected on surrounding atoms, while for 
densities n2D ~ 1 /*tm~ 2 the particle number has been 
decreased by 3% (with cooling) or 15 % (no cooling) af- 
ter a time Tg = 120 /xs. Similar to Fig. [3] we observe an 
accelerated loss of particles, resulting in a dramatic de- 



crease of the particle number for times t > Tg. The effect 
of an in-plane optical lattice is small and hardly changes 
the heating rate or particle number shown in Fig. [T4) 



2. AC-microwave field dressing 



Fig. [15] shows heating dynamics in a gas of Rydberg 
dressed atoms using the microwave dressing scheme of 
Fig. Qc) with the same parameters as in Sec. II C For 
A r = 27T x 1.0 GHz and a = 0.5 we find r c ~ 1.38 ^m. 
With uj± — 2tt x 100 kHz the hight of the barrier at the 
saddle points is ~ 70 /xK. 

In our treatment, the basic heating process of a crys- 
tal comes from its spatial rearrangement, once one of 
the atoms decays to an intermediate state |m) with zero 
dipole moment. We explain this as follows: Due to the 
lack of dipole moment, particles in the intermediate state 
|m) can travel a distance vn/T m without interacting with 
neighboring atoms, where vr is the recoil velocity and T m 
the mean lifetime of an intermediate state, (i) if vn/Y m 

is of the order of the mean particle distance in the crystal 

-1/2 



L 2D 



two particles will most probably come closer than 
the critical radius n oss and get lost. We find that this 
latter process can result in an unusual collisional evap- 
orative cooling effect, (ii) if vn/T m is smaller than the 
mean particle distance in the crystal the atom may return 
to the ground state with a finite dipole moment, resulting 
in large time-dependent fluctuations of the dipole-dipole 
interaction between neighboring atoms. The net heating 
rate of the crystal is a competition between processes (i) 
and (ii). 

We note that accidental vdW interactions may occur 
between atoms in one of the intermediate states |m), 
which may lead to additional heating. However, these 
processes are not considered here. 



Fig. 



15 



a,) shows the mean heating rate, 7g-E"kin(i = Tg) 
as a function of the atomic density for different cooling 
and lattice parameter, with -Ekin(^) the mean kinetic en- 
ergy of all N atoms. Again, each marker in Fig. [15] is an 
average over 50 runs of the simulation, where we simu- 
lated the dynamics for a time Tg = 12.8 ms. Similar as 
in the DC-electric field case of Fig. [Mj there is a strong 
dependence of the heating rate on the density and on the 
laser cooling rate /3, while the effect of an in-plane opti- 
cal lattice is weak. Additional laser cooling of the atoms 
with a rate /3 = 30 kHz again leads to a steady state of 
the mean kinetic energy (see panel (a)) and diminishes 
the particle loss rate, which is shown in Fig. |T4^b) . The 
number of remaining particles after a time Tg is shown as 
a function of the atomic density. Remarkably, more than 
70% (50%) of the particles are left with (without) laser 
cooling and densities of 122D = 0.2//im 2 after t = Tg. 
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VII. CONCLUSIONS AND OUTLOOK 

In this work we have investigated the long-time non- 
equilibrium dynamics of an ensemble of cold ground state 
atoms, weakly admixed with a Rydberg state using laser 
light. For times comparable to or larger than the effective 
lifetime of the ground state, the population of intermedi- 
ate Rydberg states following spontaneous emission signif- 
icantly affects the atomic motion in a strongly interacting 
gas, by providing a dominant heating and loss mecha- 
nism. Numerical simulations indicate that, due to the 
absence of dipole moments in the intermediate Rydberg 
states, ensembles of atoms polarized by AC-microwave 
fields exhibit a significantly different long-time dynam- 
ics, compared to atoms dressed by DC-electric fields. 

Understanding this long-time dynamics is relevant for 
the experimental demonstration of laser-dressing tech- 
niques [EI] as well as the creation of long-lived strongly 



correlated atomic phases such as self-assembled atomic 
crystals. Related to recent work on mixtures of po- 
lar molecules and Rydberg atoms [521 153], we specu- 
late that long-lived cold crystals may be used as cold 
reservoirs for achieving, e.g., sympathetic cooling of po- 
lar molecules [67] . 
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